Questions:
Does a model derived gradient between wetlands and uplands correspond to a gradient of soil development along a catena/hillslope characterized in a hydropedological framework?
If the ends of the wetland/upland classification are accurate, does the intermediate range correspond with hypothesized soil formation characteristics according to a hydropedological framework?
From Hubbard Brook website https://hubbardbrook.org/about-the-forest/:
Hydropedological units in the Hubbard Brook Experimental Forest include:
## [1] "ATTRIBUTE" "WETLAND_TY" "ACRES" "Shape_Leng" "Shape_Area"
Merge hydrography and NWI
Make terrain metrics
hbslp_3 <- SlpAsp(hbdem, w = c(3, 3), metrics = "slope", filename = "UplandWetlandGradient/data/derived_data/hbslp_3.tif",
overwrite = T)
hbslp_27 <- SlpAsp(hbdem, w = c(27, 27), metrics = "slope", filename = "UplandWetlandGradient/data/derived_data/hbslp_27.tif",
overwrite = T)
hbslp_81 <- SlpAsp(hbdem, w = c(81, 81), metrics = "slope", filename = "UplandWetlandGradient/data/derived_data/hbslp_81.tif",
overwrite = T)
hbtpi_3 <- TPI(hbdem, w = c(3, 3), shape = "rectangle", stand = "none",
na.rm = TRUE, filename = "UplandWetlandGradient/data/derived_data/hbtpi_3.tif",
overwrite = T)
hbtpi_27 <- TPI(hbdem, w = c(27, 27), shape = "rectangle", stand = "none",
na.rm = TRUE, filename = "UplandWetlandGradient/data/derived_data/hbtpi_27.tif",
overwrite = T)
hbtpi_81 <- TPI(hbdem, w = c(81, 81), shape = "rectangle", stand = "none",
na.rm = TRUE, filename = "UplandWetlandGradient/data/derived_data/hbtpi_81.tif",
overwrite = T)
hbcurv_3 <- Qfit(hbdem, w = c(3, 3), metrics = c("meanc", "profc",
"planc"), unit = "degrees", na.rm = TRUE, filename = "UplandWetlandGradient/data/derived_data/hbcurv_3.tif",
overwrite = T)
hbcurv_27 <- Qfit(hbdem, w = c(27, 27), metrics = c("meanc",
"profc", "planc"), unit = "degrees", na.rm = TRUE, filename = "UplandWetlandGradient/data/derived_data/hbcurv_27.tif",
overwrite = T)
hbcurv_81 <- Qfit(hbdem, w = c(81, 81), metrics = c("meanc",
"profc", "planc"), unit = "degrees", na.rm = TRUE, filename = "UplandWetlandGradient/data/derived_data/hbcurv_81.tif",
overwrite = T)
Slope at different scales
TPI at different scales
Curvature at different scales